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1. Introduction 

Propagation of electromagnetic or ultrasonic waves in billiards £Q , scattering of light 
in random media and transport of electrons through quantum dots [21 El share at least 
one feature in common: In all these situations one deals with an open wave-chaotic 
system studied by means of a scattering experiment, see Fig. 1 for an illustration. Here, 
we have a typical transport problem where the fundamental object of interest is the 
scattering matrix S, which relates linearly the amplitudes of incoming and outgoing 
fluxes. However, under real laboratory conditions there is a number of different 
sources which cause that a part of the flux gets irreversibly lost or dissolved in the 
environment. As a result, we encounter absorption and have to handle the S-matrix, 
which is no longer unitary. Statistics of different scattering observables in the presence 
of absorption are nowadays under intense experimental and theoretical studies. One 
should mention, in particular, experiments on energy correlations of the S'-matrix 
03 and total cross-sections [Hj , distributions of reflection E] and transmission 
[H] as well as that of the complete S matrix [H] in microwave cavities, properties of 
resonance widths |10| in such systems at room temperatures, dissipation of ultrasonic 
energy in elastodynamic billiards [TT], fluctuations in microwave networks ( see 
also references in these papers). Theoretically, statistics of reflection, delay times 
and related quantities were considered first in the strong JSj or weak ^3 absorption 
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Figure 1. A sketch of a typical experimental setup with microwave billiards. A 
flat chaotic cavity is feeded with microwaves through an attached coaxial cable 
(i.e. a scattering channel). On average, 1 — T part of the incoming flux, where 
T < 1 is the so-called transmission coefficient, is reflected back directly from the 
cable-cavity interface (port) without exciting long-lived resonances in the cavity. 
If the cavity is thin enough then only a transverse electric wave can propagate 
inside. The electric field has only a vertical component, which is uniform in 
vertical direction and distributed nontrivially in the plane. Therefore, there is 
a voltage between plates as well as a current due to the in-plane magnetic field. 
The impedance is a quantity which relates linearly the port voltage to the port 
current. Fluctuations of eigenmodes and eigenfrequencies result in fluctuations of 
the impedance or S-matrix, as the driving frequency or port position is changed. 



limits at perfect coupling, and very recently at arbitrary absorption and coupling 

tuminiEiEU]. 

Another insight to the same problem comes by considering it not from the 
"outside" , but rather from "inside" . Then the impedance relating linearly a voltage 
to a current turns out to be the prime object of interest |2l)l 121) . see Fig. 1. It turns 
out that after proper taking into account of the wave nature of the current |221 12~5| , 
the cavity impedance becomes an electromagnetic analogue of Wigner's reaction (R) 
matrix of the scattering theory. This can be understood qualitatively through the 
well-known equivalence of the two-dimensional Maxwell equations to the Schrodinger 
equation, the role of the wave function being played by the field (the voltage in our 
case). Then the definition of the impedance becomes formally similar to the definition 
of the i?-matrix (which relates linearly the normal derivative of the wave function to 
the wave function itself on the boundary). The impedance is, therefore, related to 
the local Green function of the closed cavity and fluctuates strongly due to chaotic 
internal dynamics. 

The imaginary part of the local Green function (which is proportional to the real 
part of the impedance) is known in the context of mesoscopics as the local density of 
states and has a long story of study, see [21] for a recent review. Actually, a closely 
related quantity in the context of spectra of complex atoms and molecules has the 
meaning of the total cross-section of indirect photoabsorption, see e.g. |2S]- As to the 
real part, it seems to have no direct physical meaning in mesoscopics while it has the 
meaning of reactance in electromagnetics, where both real and imaginary parts are 
experimentally studied. Very recently an approach |2fil I27j has been developed by us 
which allows one to study the (joint) distribution function of these real and imaginary 
parts at arbitrary absorption and to relate it to the reflection distribution, thus linking 
somewhat complementary experiments jS] and |20j together. 

Due to a strong resonance energy dependence the impedance and 5-matrix as 
well as any scattering observable exhibit strong fluctuations over a smooth regular 
background as the scattering energy (or another external parameter) is varied. These 
two variations occurring on different energy scales are usually decomposed into a mean 
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and a fluctuating part by means of the spectral or (assumed to be equivalent) ensemble 
average (•■■}. In this paper we consider statistics as determined by a two-point 
correlation function of the fluctuating parts (also called a "connected" correlator): 
{AB) conn — (AB) — (A) (B). We restrict ourselves below to the cases of preserved and 
broken time-reversal symmetry (TRS). 



2. Scattering, RMT and absorption 

The resonance energy dependence of observables becomes explicit in the well-known 
Hamiltonian approach to quantum scattering, which was developed first in the context 
of nuclear physics [2E1 HID EH an d adopted later for the needs of mesoscopic physics, 
see e.g. This framework is adequate to take finite absorption into account 

as well. We have the following relation between the resonance part of the scattering 
matrix and Wigncr's reaction matrix: 

S{E)= \V\K{E) ' K{E) = \V\E-H)-^V. (1) 

The Hamiltonian H of the closed system gives rise to N levels (eigenfrequencies) which 
are coupled to M continuum channels via the NxM matrix V of coupling amplitudes 
V£ (n — 1 • ■ ■ N, c = 1 . . . M). Performing for S a Taylor series expansion in K and 
regrouping the terms, one comes to another well-known expression for the 5* matrix 

S(E) = l-W \ 1 V, H eS =H-iVVt. (2) 

in terms of the effective Hamiltonian Ti c s of the open system, which is non-Hermitian 
contrary to the Hermitian H. The factorized structure of the anti-Hermitian part 
ensures the unitarity of S(E) at real values of E. In a resonance approximation of the 
energy-independent amplitudes the complex eigenvalues £ n — E n — ^T n of H s are 
the only singularities of the S matrix in the complex energy plane. As required by 
causality , they are located in the lower half plane and correspond to the long-lived 
resonance states, with energies E n and escape widths r„ > 0, which are formed on 
the intermediate stage of a scattering process. 

To mimic chaotic nature of the intrinsic motion we adopt, as usual, the random 
matrix theory (RMT) 031 IH3 and replace the actual Hamiltonian with a random 
Hermitian matrix H. It turns out that spectral fluctuations possess a large degree of 
universality in the limit N—>oo: being expressed ("unfolding") in units of the mean 
level spacing A they become independent of microscopic details (i.e. a particular 
form of the distribution of H) and get uniformly distributed throughout the whole 
spectrum |35) . That amounts usually to considering local fluctuations at the center 
of the spectrum (E = 0) and to restricting ourselves to the simplest case of Gaussian 
ensembles. On has the Gaussian Orthogonal Ensemble (GOE, the Dyson's symmetry 
index f3 = 1 and H symmetric) for chaotic systems with preserved time-reversal 
symmetry (TRS) and the Gaussian Unitary Ensemble (GUE, (3 — 2 and H Hermitian) 
for those with fully broken TRS. For similar reasons the approach is independent of 
particular statistical assumptions on coupling amplitudes V® (as long as M <C -/V 
j.'itil I37p. which may be chosen as fixed [2H| or random variables. They enter 
final expressions only by means of M transmission coefficients (also so-called sticking 
probabilities) 

r -i iggia- ^ K _?WT (3) 

lc ~ l lb 1 " (1 + Kc ) 2 ' Kc " 2iVA ■ {3) 
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where S cc stands for the average (or "optical") S matrix. They are assumed to be 
input parameters of the theory. T c <C 1 or T c = 1 corresponds to an almost closed or 
perfectly open channel "c" , respectively. 

Absorption is usually seen as a dissipation process, which evolves exponentially 
in time. Strictly speaking, different spectral components of the field have different 
dissipation rates. However, this rather weak energy dependence can easily be neglected 
as long as local fluctuations on much finer energy scale ~ A are considered. As a 
result, all the resonances acquire additionally to their escape widths one and the same 
absorption width T > 0. The dimensionless parameter 7 = 27iT/A characterizes then 
the absorption strength, with 7 C 1 or 7 > 1 corresponding to the weak or strong 
absorption limit, respectively. (Microscopically, it can be modelled by means of a 
huge number of weakly open parasitic channels EHj or by additional coupling to 
very complicated background with almost continuous spectrum \T^, see also 

Treating T phenomenologically, one sees that such a uniform absorption can 
equivalently be taken into account by a purely imaginary shift of the scattering energy 
E — > E + ±T = E 7 , so that the ^-matrix S~,(E) = S(E^) becomes subunitary. The 
reflection matrix SJ^S^ < 1 provides then a natural measure of the mismatch between 
incoming and outgoing fluxes |141 115| . At last but not least, the matrix Z = iK(Ej) 
has the meaning of the normalized cavity impedance in such a setting, see for 
further details. 



3. Correlation functions 

3.1. Impedance 

Let us consider first the simplest case of the impedance when the problem can be fully 
reduced to that of spectral correlations determined by the two-point cluster function 
Y 3 ,p((d) = S(uj) - A 2 (p(Ei)p{E 2 )) conn , where u = (E 2 - E^/A and p(E) being the 
spectral density, ft is easy to find (Z ab ^ = K a S ab for the mean impedance at E = 0. 
To calculate the energy correlation function 

C a z bcd (ou) = {Z ab *{E 1 )Z cd {E 2 )) conn (4) 

it is instructive to write Z ab (E) — | J2n e-e +Tr/2 m tne e ig en basis of the closed 
system. The rotation that diagonalizes the random H transforms the (fixed) coupling 
amplitudes V£ to gaussian distributed random coupling amplitudes v® with the zero 
mean and the second moment = (2K a A/7r)i5 ab (5 rim . In such a representation 

10} acquires the following form: 

C abcd {u) = y\ (vy b n *v%v±) ( ; =- ) (5) 

Z y ' f^4\ " " m m/ \ El -E n -^TE 2 -E m + ^T / yJ 

n,m \ * ^ I conn 

so that averaging over coupling amplitudes (i.e. eigenfunctions) and that over the 
spectrum can be done independently. The gaussian statistics of v results in 

1 / 7r \2 / a &* c* d \ cab zed , / rac cbd , r cad ?bc\ r fr>\ 

i(^) {v n v n v m v m ) = K a K c b 6 + K a Hb{0 +difjd 6 jd nm (6) 

where 5ip term accounts for the presence of TRS, when all wJJ are real and Z is sym- 
metric. It is useful then to represent the spectral correlator in the form of the Fourier 
integral dt, / °° ^e-^+^^e 1 ^-^ ^i(E ntl -E m t 2 )^^ Due 
to the uniformity of local fluctuations in the bulk of the spectrum, one can integrate 
additionally over the position E of the mean energy: J -^-e 1 ^* 2- * 1 ) = -^■(5( * 2 t ~ tl ), 
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Figure 2. The impedance (a) and S-matrix (b) enhancement factors for chaotic 
systems with preserved (/3 = 1) or broken (/3 = 2) time-reversal symmetry as 
functions of the absorption strength 7. The case (b) corresponds to the many- 
channel limit with isolated Q2 C '^ C = son< ^ lines) or overlapping resonances 
T c = 5, dashed lines). 



where tn = 2tt/ A is the Heisenberg time. From the known RMT spectral fluctuations 
one also has (1 — JV) (e 1 ^" - ^'')^, = &2,/3(i/te) for n^m, where b2,p(r) is the 
spectral form factor defined through the Fourier transform of Y2„s(w) |341 135j : 

6a,/9=i(r) = [1 - 2t + r log(l + 2t)]0(1 - r) + [t log(^) - 1]6(t - 1) (7a) 
&2,/3= 2 (t) = (1 - r)6(l - r) (76) 

at t > and b2,p{—r) = o 2 ,p{t), so that Y 2 ^{uj) — /^ o dTe 27rl " r 6 2 ,/3(T). Combining 
all these results together and measuring the time in units of the Heisenberg time 
(t = t/tn), we arrive finally at 



Cf ca {w) = / dre ZmulT Cf ccl {T) (8a) 
Jo 

Cf cd {r) = 4 e-^ [k k c (1 - 6 2:/3 (r))5 ah ( 5 C!i + K a K b (6 ac 6 bd + 5 10 5 ad 8 bc )\ . (86) 

Similar in spirit calculations were done earlier a context of reverberation in complex 
structures in ^2 an d in a context of chaotic photodissociation in |421 BS] . 

The form factor (|8"&"|) is simply related to that of K matrix elements at 
zero absorption as Cf bcd (r) = e _7T C^- fccd (r). Such a relationship between the 
corresponding form factors with and without absorption is generally valid for any 
correlation function which may be reduced to the two-point correlator of resolvents 
(see jS] and below, e.g., for the case of the S matrix). This can be easily understood 
as the result of the analytic continuation 2t:lj — > 2t:lj + i"f of the energy difference ui 
when absorption is switched on (see the previous section). 

The obtained expressions describe a decorrelation process of the Z matrix 
elements as the energy difference grows, generally, Cz(u— >oo) — ► 0. At co = 0, 18 all 



provides us with impedance variances Cf ab (Q) = var(Z afc ) = {\Z ab \ 2 ) - \{Z ab )\ 2 , 
which were recently studied in 03] ( see also UM ) • I n analogy with the so-called elastic 
enhancement factor considered frequently in nuclear physics |46| . one can define the 
following ratio of variances in reflection (a=6) to that in transmission (a^=b) : 
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where the second equality follows easily from l|86(l (note that the coupling constants 
K a> fc are mutually cancelled here). Making use of & 2 ,/3(oo)=0 and 6 2j/ g(0)=l, one can 
readily find Wz,p in the limiting cases of weak or strong absorption as: 

at 7 < 1 

10 

at 7 > 1 

Wz,p decays monotonically as absorption grows, see Fig. la. In the case of unitary 
symmetry, (|76|l and 10 yield explicitly Wz,2 = — e~ 7 ) in agreement with 03]. It 

is hardly possible to get a simple explicit expression at finite 7 in the case of orthogonal 
symmetry However, a reasonable approximation can be found if one notices that 
the integration in is determined mainly by the region s < 1, so that one can 
approximate &2,i(s) ~ (1 — 2s + 2s 2 )0(l — s) through its Taylor expansion. Performing 
the integration, one arrives at Wz,i ~ 3 — 7~ 2 [(4 + 7 2 )(1 — e / ) — 27(1 + e *)], which 
turns out to be a good approximation to the exact answer at moderate absorption 
(deviations are seen numerically only at 7 ~ 1). 



3.2. S -matrix elements 

The energy correlation function of the scattering matrix elements 

/>OC 

Cg cd (« J ) = (S?'(E 1 )S?(E 2 )) conu = / dre^Cf cd {r) (11) 

Jo 

is a much more complicated object for an analytical treatment as (0J. The reason 
becomes clearer if one considers again the pole representation of the S matrix which 
follows from ©: S ab (E) = S ab - iYjW^w b J(E - £„). Due to a unitarity constraint 
imposed on S (at real E), the residues and complex energies get mutually correlated 
|3(J| with a generally unknown joint distribution. The separation like JSJ into a 
"coupling" and "spectral" average in no longer possible and can be done only by 
involving some approximations |47| . The powerful supersymmetry method |291 148| 
turns out to be an appropriate technique to perform the statistical average in this case. 
In their seminal paper |29j . Verbaarschot, Weidenmiiller and Zirnbauer performed the 
exact calculation of (|11J) at arbitrary transmission coefficients (and zero absorption) 
in the case of orthogonal symmetry. This finding was later adopted [S] to include 
absorption. The corresponding exact result for unitary symmetry is still lacking in 
the literature (see, however, 02] concerning the iS-matrix variance in the GOE-GUE 
crossover at perfect coupling) and will be presented below. 

The calculation proceeds alone the same line as in |29| , we indicate only essential 
differences. As usual, the representation of resolvents and thus (|llfl in the form of 
gaussian integrals over auxiliary "supervectors" consisting of both commuting and 
anticommuting (Grassmann) variables allows one to perform statistical averaging 
exactly. In the limit N — > 00, the rest integration over the auxiliary field can be 
done in the saddle-point approximation. The final expression for both the correlator 
and its form factor (|11|) can be equally represented as follows: 

gated = S ab 5 cd TaTc ^l _ Ta)(1 _ + ( S ac § bd + Sw S ad S bc )T a T b P ab . (12) 

Here, the Sip term accounts trivially for the symmetry property S ab — S ba in 
the presence of TRS. J ac and P a b defined below are some functions (of the energy 
difference to or the time r), which depend also on TRS, coupling and absorption 
but already in a nontrivial way. As a result, the elastic enhancement factor 
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Ws,p = \J va.i{S aa )v'&v(S hb ) /vax{S ah ) is generally a complicated function of all these 
parameters, in contrast to 0. In the particular case of perfect coupling, all T c = 1, 
one has obviously from Ijl2|l that Ws,p = 1 + Sip at any absorption strength. 

The saddle-point integration turns out to have a nontrivial saddle-point manifold 
[J5J over which one needs to integrate exactly. This task can be accomplished by 
making use of the "angular" parametrization [2H| of the manifold in terms of the ^ x -| 
supermatrices ti2 and £21- We consider first real u> (no absorption). Then the functions 
J ac and P a b have in the energy domain the following representation: 







1G 

(3 



t-\iti 



T a tl2t2l 



^21 t-i 



J ac (uj) = '- ( str - \ k str - k J&> ) (13a) 



1 + T c t2ltl2 



1G 



(i + ti 2 t 21 )V 2 ^ (l + t 2 iti 2 ) 1/2 ,j 



p Qb ( W ) = ^ (str t 21 \',^v ; ^ \ ,:v ; fc j~m J (V.\U) 



1 + 7a*12*21 1 + T b t 2 lti2 , J( 

that is completely in a parallel with (2H| (the diagonal matrix fc = 1 (—1) in the 
subspace of commuting (anticommuting) variables). This result has the form of an 
expectation value ((• • = fd/j,(p) (• • -)e lCl3 ^' in the field theory (nonlinear "zero- 
dimensional" supersymmetric er-model) characterized by the Lagrangian Cp{uj) = 

Pttlu str (t 12*21)- The so-called channel factor = nfli S( f e t(l + Tcii2*2i)~^ 2 

accounts for system openness. We refer the reader to [23 ISHI for a definition of the 
supertrace and superdeterminant as well as for a general discussion of the superalgebra. 
An explicit parametrization of matrices ii 2 , i 2 i and the integration measure d/Lt(^) over 
them depend on the symmetry case considered; it can be found in |29| for (3 = 1 and 
in |^ for (3 — 2. Essential is that the final expressions are determined only by 
real "eigenvalues" /io and /ii. 2 of the angular matrices. Finally, one can cast resulting 
expressions as follows 

JM = ((j^ + ife + i*;) (1*7 + rfe + i*r)^) p ( 14 «) 

P ■ (,,,\ — It Ml+M I M2(l+M2) I ^o(l-Mo) \T,r\ (~\AK\ 

r a b\U) - \l(i +TaMl )(l+T 6Ml ) + (l+T aM2 )(l+T 6M2 ) + (l-T^Kl-T^oji^/,, 

with J r M = IIc[ (i+t 1 ^"i" Z )''(i+t * n * ne ^ = ^ case °^ or fhog ona l symmetry 

and 



JaeH = ( (T+ffer + rfe) (rfc + t^)^m ) m (15a) 

with JF^/ = Yl c i+y c ^° m f ne /3 = 2 case of unitary symmetry. Here, the corresponding 
integration ((• • is to be understood explicitly for these two respective cases as 

u d/ii X dM2 io dM0 id + + ^m»o + M!) 2 (^ + ^ ( ■ (16) 

and 

fl 

-2 i27rw(/Ji+/i ) 



d/xx / d// (m + Mo)" 2 e ^^i+«J (...). (17) 
Jo 

In the important particular case of the single open channel (elastic scattering), the 
general expression for the /3 = 2 case simplifies further to 

(S*(Ei)S(E 2 )) coan = T 2 /°V f -^- ^ilruV 1 ^ 1+tl ° ] • (18) 
Jo Jo Mi + Mo U + -* Mi J 
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Finally, putting above ui — > to + vy /2ir accounts for the finite absorption strength 7. 

To consider ()1 2t in the time domain, i.e. the form factor Cg bcd (r), we notice that 
the variable r = \{^\ + ^2 + 2/io) for /3 = 1 or r = + ^ for P — 2 plays the role of 
the dimensionless time. The corresponding expressions for P a b{T~) and J ac (r) can be 
investigated using the methods developed in 0B| 03 E3 . For orthogonal symmetry 
it was done in j^j, where the overall decaying factor e _7T due to absorption was also 
confirmed by comparison to the experimental result for the form factor measured in 
microwave cavities. It is useful for the qualitative description to note that P a b{T~) 
and 2J ac (r) are quite similar to the "norm leakage" decay function [SJ and the form 
factor of the Wigner's time delays %iT\, respectively (they would coincide exactly 
at 7 = 0, if we put T a ^ iC = appearing explicitly in denominators of i|13a[) and 
thereafter). Then one can follow analysis performed in these papers, see also |53j . 
to find qualitatively P ab {r) ~ e^l + f^rT^l + f^r)- 1 n c (l + ^T cT y^ 2 and 
Jab(i~) ~ [l-6 2 , / 3(r)]P Qb (r). One has P afc (r) rj e~ 7T and J afc (r) « (2T//3)e" 7r as exact 
asymptotic at small times 07], they both being ~ e -7T T -M0/2-2 a ^ i ar g e times. 

Such a power law is characteristic for open systems [21 ESI EH • Physically, it 
results from width fluctuations, which diminish as the number M of open channels 
grows |22 I n the limiting case M — > 00 and T c — ► 0, all the resonances acquire just 
the same escape width (in units of t^ 1 ) 5Z C T C , which is often called the Weisskopf's 
width [nS], so that the total width is 7t = ^2 C T C + 7- Then further simplifications 
occur: P a b{T) = e~ lTT and J a b(r) = [1 — 62, /^t)]^^ 7 ", that results finally in 

Cf cd {u) = {daCSbd + +S ab S cd T a T c / 00 dr[l-6 2 M^)]^ (7r ~ 2 ^ )T -(19) 

7T - 27Tio; J Q 

For the case of (3 = 1 this result (at zero absorption) was obtained earlier by 
Verbaarschot 021 ■ In the limit considered, expression i|19[) is very similar to (|8af) . 
II8&L so that the enhancement factor Ws.p is given by the same (0 where 7 is to 
be substituted with jt, see Fig. 1(b) for an illustration. At 7y 3> 1 (large resonance 
overlapping or strong absorption, or both) the dominating term in Ijl9(l is the first one, 
which is known as the Hauser-Feshbash relation jHE], see j^ZI EH1 EH] for discussion. 
Then Ws,p = 2/(3 = Wz,p that can be also understood as the consequence of the 
gaussian statistics of S (as well as of Z) in the limit of strong absorption 59 . 

4. Conclusions 

For open wave chaotic systems with preserved or broken TRS we have calculated 
exactly the energy correlation function of impedance matrix elements at arbitrary 
absorption and coupling. This function is found to be related to the two-level cluster 
function, or to its form factor in the time domain. The overall exponential decay due 
to uniform absorption is shown to be the generic feature of any correlation function 
reduced to a two-point spectral (resolvent) correlator, that follows simply from analytic 
properties of the latter in the complex energy plane. The elastic enhancement factor 
defined though the ratio of variances in reflection to that in transmission diminishes 
gradually from the value 1 + 2/(3 at weak absorption to 2/(3 at strong absorption. 

The similar exact calculation for S'-matrix elements has been performed in the 
case of broken TRS, thus completing the well-known result [211 of preserved TRS. The 
corresponding enhancement factor never reaches the maximum value 1 + 2//3 at any 
finite resonance overlapping. It attains the value 2/(3 in the limit of strong absorption 
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(independent of coupling) or at perfect coupling (independent of absorption). 
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